Pulse calibration and non-adiabatic control of solid-state artificial atoms 
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Transitions in an artificial atom, driven non-adiabatically through an energy-level avoided cross- 
ing, can be controlled by carefully engineering the driving protocol. We have driven a superconduct- 
ing persistent-current qubit with a large-amplitude, radio- frequency field. By applying a bi-harmonic 
waveform generated by a digital source, we demonstrate a mapping between the amplitude and phase 
of the harmonics produced at the source and those received by the device. This allows us to image 
the actual waveform at the device. This information is used to engineer a desired time dependence, 
as confirmed by detailed comparison with simulation. 

PACS numbers: 03.67.Lx, 32.80.Qk, 78.70.Gq, 85.25.Cp, 85.25.Dq 



Due to the strong coupling of solid-state artificial 
atoms []], Q to electromagnetic fields, a wide range of 
pulsing techniques can be used to manipulate and con- 
trol such systems. A common approach, first developed 
for natural atomic systems, involves low- frequency Rabi 
oscillations driven by low-amplitude microwave radiation 
3, 0, H, Ej, combined with envelope-shaping techniques 
3, |8[. However, strong driving affords many more possi- 
bilities for quantum control. Fast dc pulses can be used 
to bring the system non-adiabatically into the vicinity 
of an avoided level-crossing, where coherent oscillations 
result from Larmor-type precession [9I. [ToL Hl|. 



Furthermore, quantum coherence was recently shown 
to persist under lar ge-a mplitude harmonic driving [l2|, 



IS UJ, Ha llfl HI UM' In this regime, by carefully 



engineering the driving protocol, arbitrary rotations of 
a qubit's quantum state on the Bloch sphere could be 
performed. Such protocols may lead to much faster 
quantum-logic gates than those achieved by using Rabi- 
transition based techniques. To make such an approach 
feasible, one must be able to apply external fields of ar- 
bitrary time dependence to a quantum device. With this 
ability, the techniques used for pulsed NMR can be ex- 
tended to achieve even higher-fidelity quantum control 

In addition to the challenge of designing optimized 
pulses, accurately delivering them to a device in a cryo- 
genic environment is a difficult problem in its own right, 
especially at radio and microwave frequencies. In par- 
ticular with transient pulses @, E3, El , it is hard to de- 
termine the exact pulse shape at the device. Although 
digital waveform generators offer the possibility to cre- 
ate control pulses with essentially arbitrary time depen- 
dence, the signal that reaches the device may be strongly 
distorted due to the impedance mismatch and frequency 
dispersion that can occur in the long coaxial cables lead- 
ing from the generator to the device. Thus, in order to 



achieve high-fidelity quantum control, it is important to 
learn not only how to design optimized pulses with arbi- 
trary shape in the computer, but also how to apply them 
faithfully to the device under study. 

In this Letter we present methods to image the ex- 
tremal points of a periodic signal and the actual wave- 
form of a single pulse, as received by the device. Ap- 
plied to digitally generated, radio-frequency bi-harmonic 
pulses, Eq. (pQ), this approach allows us to determine the 
amplitude ratio and phase between harmonics, as illus- 
trated in Figs [1] and 03 Using this information, we create 
waveforms that are calibrated to accurately produce the 
desired time dependence at the device. In particular, the 
ability to engineer waveforms is used to control the time 
spent near an avoided crossing. The effect of the lat- 
ter on quantum evolution is illustrated and confirmed by 
detailed comparison with simulations. 

In our experiments, we use a niobium superconduct- 
ing persistent-current (PC) qubit, see Fig. QJa). When 
the magnetic flux / piercing the qubit loop is nearly 
$o/2, where $0 = h/2e is the flux quantum, the qubit's 
potential energy exhibits a double- well profile. The 
system then supports a set of discrete, diabatic states 
L ), I g, R); (p, q) = 0, 1, . . .}, localized in the left and 
right wells, and associated with opposing persistent cur- 
rents. Their energies vary linearly with the magnetic flux 
and exhibit avoided level crossings, see Fig. QJb). 

We drive the system longitudinally with an applied 
flux 5f(t) = f(t) - $ /2 = 5/dc + AfO) consisting of 
a static and a time-dependent part. This driving cre- 
ates large-amplitude excursions through the energy-level 
diagram from the point of the static flux bias % c , as in- 
dicated by the arrows in Fig. QJb). The quantum evolu- 
tion is primarily adiabatic, except when the field sweeps 
through the energy-level avoided crossings A p?g where 
Landau-Zener (LZ) transitions split the state into coher- 
ent superpositions of | p, L ) and \q,R). When the excur- 
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FIG. 1: (color online) (a) Schematic experimental set-up. The 
device was placed at 20 mK (see [13] for details), (b) Energy 
diagram determined by amplitude spectroscopy [16[ . The en- 
ergy gaps are small on the scale of the figure (Ao, q /h = 0.013, 
0.090, 0.40, and 2.2 GHz for q = 0, . . . , 3, resp.). Below: 
flux excursion Sf(t) during the short, bi-harmonic pulse that 
probes the time-dependent dynamics in Fig. [3] (c-d) Exam- 
ples of bi-harmonic waveforms JT]). values are indicated by 
dashed, vertical lines; blue/red denote maxima/minima. The 
lower left waveform corresponds to the one in Fig. EJb). 



sion extends over multiple levels, the interplay between 
LZ transitions at different crossings and intrawell relax- 
ation (\p',L) — > \p,L) and \q f ,R) — > \q,R)) gives rise 
to "spectroscopy diamonds" in the saturated population 
[l6[ as shown in Fig. [2(a). Because the LZ-transition 
probability Plz = 1 — exp(— irAp q /2hv) is exponentially 
sensitive to the energy sweep rate v through the avoided 
crossing, the magnitude of population transfer contains 
useful information about the driving signal in the vicinity 
of the avoided crossing. The bright diamond edges occur 
for combinations of static flux detuning Sf^c and rf am- 
plitude where an extremum of the driving waveform just 
reaches an avoided crossing. At these parameter values, 
the system spends the most time in the vicinity of the 
avoided crossing, thus allowing for maximum population 
transfer. As described below, this phenomenon allows us 
to use the bright bands due to transitions at a particular 
avoided crossing to image the extrema of an rf pulse. 

In our experiment, we synthesized the control pulses 
digitally using the Tektronix AWG5014 arbitrary wave- 
form generator with effectively 250 MHz analog output 
bandwidth. The pulses were launched through a coaxial 
cable into an on-chip coplanar waveguide terminated in 
an "antenna" (loop of wire). The frequency-dependent 
impedance of this antenna, as well as the resistive coaxial 



cable, makes our microwave line differ from an ideal 50- 
Q system. Moreover, its a priori frequency-dependent, 
inductive coupling to the qubit calls for amplitude and 
phase calibration at the device for each frequency com- 
ponent of the pulse. 

To demonstrate our method of pulse imaging and cal- 
ibration, we drive the qubit with the bi-harmonic signal 



/ rf (t) = A 1 cos(cjt+#) + A 2 cos(2ujt). 



(1) 



We classify the waveforms, described by Eq. (pQ), in terms 
of the relative amplitude ratio and phase between har- 
monic components using the notation {A2/A1, 0}. The 
waveforms can have either two or four extrema per cycle, 
depending on the values A2/A1 and see Figs QJ c-d). 

The cross-over between the domains of two and four 
extrema per cycle occurs at A2/A1 = 1/4. At this am- 
plitude ratio, the phase = tt is particularly interesting: 
this waveform has one parabolic (quadratic) and one flat 
(quartic) extremum, the flatness of which arises from can- 
celation of the quadratic terms in the Taylor expansions 
of the two harmonic components. The presence of the 
flat extremum relies on the delicate balance of amplitudes 
and phases; hence we exploit the {1/4, 7r}-waveform to 
calibrate the relative phase difference between the two 
frequency components at the device as shown in Fig. O 

Each experiment uses the following pulse sequence, see 
Fig. QJa). A static flux is applied using a superconduct- 
ing coil, and maintained fixed at some value 5fd c - A har- 
monic, 11-MHz, 3-/is cooling pulse initializes the qubit to 
its ground state [20[. After a short delay (~20ns), an rf- 
driving pulse of duration At and amplitude / r f is applied. 
The mixer in Fig. Ufa) is used to terminate the pulse after 
a time At in the experiments shown in Figs 09(b,d), but 
is not used in the experiments of Fig. [2l After another 
short delay, a sample-and-hold current pulse is applied 
to the hysteretic SQUID magnetometer [2l|, which mea- 
sures the qubit in the basis of the diabatic states, and 
is set to switch only if the qubit is in state \L). The re- 
sulting voltage signal is amplified at room temperature 
and recorded with a threshold detector. This sequence is 
repeated a few thousand times at a rate of 10 kHz. 

First, we independently calibrate the amplitudes of the 
two components 1 and 2 by finding their points of optimal 
cooling [20[ in a single-frequency spectroscopy diamond, 
see Fig. EJa). These points occur when the flux excur- 
sion precisely reaches both opposite side-avoided cross- 
ings Ao,i and Ai 5 o at zero dc bias. We denote these am- 
plitudes by A\ 2 , and define A*=A\+A* 2 and A = A X +A 2 . 
By adding the two harmonics together, fixing the ampli- 
tude at a value within the range < A/ A* < 0.5 (see 
Fig-Efb), dashed line), and scanning over flux detuning 
Sfdc an d phase g at the generator, we obtain Fig. [Sfc). 
As mentioned above, high intensity occurs when an ex- 
tremum of the waveform just reaches an avoided crossing. 
Thus the value of g at the cusps in Fig. [21(c) corresponds 
to the desired value of = tt at the device [cf. Fig. [TJd)]. 
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FIG. 2: (color online) Calibration of rf amplitude and phase using long pulses, (a) "Spectroscopic diamonds" from single- 
harmonic driving. The points A^ 2 of optimal cooling (arrows) are used for amplitude calibration at each frequency, (b) 
Diamonds when driving with the bi-harmonic waveform {1/4, tt} at uj/2tv — 45 MHz, see Eq. (pQ). The arrows point out the 
diamond edges corresponding to the quartic extremum reaching the Ao, q crossings and the quadratic the A Pj o- (c) Measurement 
along the dashed line in (b) for varying g . The quartic extremum — n is obtained at a cusp (solid white line), cf. Fig. [TJd). 



Applying the calibrated, composite {1/4, 7r}-waveform 
to the device, we obtained the bi-harmonic spectroscopy 
diamonds in Fig. EJb). The diamonds are skewed com- 
pared to (a) due to the asymmetry of the signal f T f(t) 
about zero; a quadratic extremum has a larger amplitude 
than a quartic. Additionally, the data exhibit greatly en- 
hanced population transfer at the diamond edges arising 
from the flat extremum relative to that at the edges aris- 
ing from the parabolic extremum; the same is seen at the 
cusps in Fig.[2jc). As discussed above, the reason is that 
the flat maximum makes the system spend more time at 
the avoided crossing than does the parabolic maximum. 

The data in Fig. [2] reflect the averaged response of the 
qubit after a 3-/is pulse of base frequency o;/27r = 45 MHz 
(135 cycles). However, on short time scales, the pulse re- 
quires a few cycles to reach its steady- state amplitude, at- 
tributable to the build-up time of the standing- wave volt- 
age at the imperfectly terminated end of the microwave 
line. In order to adjust the pulse shape for short pulses, 
we investigate the brief, temporal Larmor-type oscilla- 
tions induced by a single, non-adiabatic Ao, 2-passage, 
and the build-up of Stiickelberg oscillations occurring 
during the first return trip through that crossing. 

When the system traverses the crossing, the LZ process 
creates a quantum superposition of the states associated 
with different wells. Upon return, these two components 
interfere with the relative phase 13, 3, 13] 



A0 
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dt Ae(t), 



(2) 



where Ae(t) is the instantaneous diabatic energy- level 
separation. This leads to Stiickelberg-interference fringes 
in the occupation probabilities, provided that the time of 
the excursion, typically a fraction of the drivi ng p eriod, 
is smaller than the relevant decoherence times [17j . 

We apply a nominally flat-top pulse with an amplitude 
that reaches through Ao,2- At a time At, we abruptly 
terminate the pulse, and so obtain a snapshot recording 
of the population, see Fig. 09(b). Experimentally this is 



done in the following way. A microwave cooling pulse ini- 
tializes the system to the ground state |0, L) at 5fa c < 0. 
Then, a bi-harmonic pulse (pQ) starts the flux excursion 
towards the left, drawing the system deep into its ground 
state, as shown by the arrows in Fig. [2(b) (in this way we 
minimize the effects of unavoidable transients associated 
with the pulse turn-on). Upon return towards the right, 
the system rapidly traverses the avoided crossings A , o 
and A ,i before reaching A ,2 [22[. We interrupt the 
excursion at the time At by mixing the pulse with the 
falling edge of a square pulse, thus returning the flux de- 
tuning to the value 5fd c - The population can then be read 
out well before inter- well relaxation occurs. Through this 
pulsing scheme, we can increment the pulse length At in 
steps smaller than the AWG's sampling time. 

For pulse lengths At long enough for the system to 
complete a round trip through Ao,2, the final popula- 
tions in the left and right wells are stationary, i.e. inde- 
pendent of At. The spacings of the resulting horizontal 
Stiickelberg interference fringes in Fig. [31(c) are deter- 
mined by the phase ([2]) gained over the course of the 
excursion beyond Ao,2; as described below, these fringes 
thus constitute a signature of the part of the waveform 
extending beyond the avoided crossing. 

The quantization condition for interference, A</>i2 = 
2nN gives the static flux detuning Sf^ for the A^ th 
node of constructive interference. For generic pulses with 
extrema that can be fit to parabolae, Eq. ([2|) leads to an 
TV 2 / 3 power law for the values of dc-flux detuning at the 
interference nodes [16[. However, for the {1/4, tt} flat- 
top waveform that we seek, the quartic extremum gives 
rise to an TV 4 / 5 power law, which is the signature we use 
to determine when the desired pulse has been attained. 
In Figs09(f-g) we demonstrate the iV 4 / 5 power law for a 
pulse that is calibrated by comparison to a simulation, in 
which we numerically integrate the Schrodinger equation 
for a two- level system with time-dependent driving (pQ). 
The initial state is taken to be the ground state of the 
system at detuning Sfdc- For further details, see pH ]. 
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FIG. 3: (color online) Time- dependent response, (a-b) Population in state \L) at the final time At of a variable- length, flat-top 
pulse {1/4, 7r}, reaching through An, 2 as in Fig. [TJb). The amplitude A is indicated by the rightmost dotted line in Fig. [2{b). 
The edges of the curved fringes provide a direct image of the pulse (black, dotted line), obtained experimentally by comparing 
(a) and (b). (c) Population in state \L) after a return trip (At = 22ns) through An, 2. (d-e) Like (b-c) but with a waveform 
calibrated using the long pulses in Fig. [2] thus failing to give the correct waveform, (f) The lower (red) and upper (blue) sets 
of data show the spacings of the Stiickelberg-interference fringes in (c) and (e), respectively, with fits to the N 4 ^ 5 power law 
(fringes N = 4 — 35 and III — XXI fitted), (g) The residuals indicate that the fringes in (e) [blue x] do not fit well, while those 
in (c) [red +] fit better; this confirms that the calibrated pulse better approximates the {1/4, ty} waveform. 



Overall there is a good agreement between simulation 
and experimental data. We attribute the discrepancy 
around At = 18 ns, Sf^c — -5m$o in Fig. [31(b) to non- 
ideal pulse turn-off, and the loss of fringe contrast to 
intrawell relaxation, digital noise, and phase jitter [16[. 

We note that a rectangular-pulse waveform would ide- 
ally lead to a very different fringe pattern compared to 
the one in Fig. [3l showing maximal contrast when the 
pulse height matches the level-crossing position in en- 
ergy, and symmetrically decreasing at higher and lower 
pulse amplitudes. Interestingly, in the pioneering work 
of Nakamura et al. [§], where such pulses were used to 
drive a qubit, the observed fringe pattern had a strong 
asymmetry, overall resembling the data in Fig. [3l This 
behavior can be attributed to a smooth turn-on and turn- 
off of the pulses used in that work. 

In this work, we present an approach that in principle 
can be used to design an arbitrary waveform at the de- 
vice, and use it to drive a superconducting qubit. This 
technique, which is demonstrated with biharmonic driv- 
ing, relies on calibrating the different Fourier harmonics 
comprising the waveform. Having performed such a cali- 
bration, we are able to control the time spent by the qubit 
near an avoided crossing. The details of qubit response 
are accounted for by comparison with a simulation. 

We gratefully acknowledge discussions with and sup- 
port from T. P. Orlando. We thank Y. Nakamura, F. 
Nori and T. Yamamoto for helpful discussions; G. Fitch, 
V. Bolkhovsky, P. Murphy, E. Macedo, D. Baker, and 
T. Weir at MIT Lincoln Laboratory for technical assis- 



tance. This work was sponsored by the U.S. Government 
and the W. M. Keck Foundation Center for Extreme 
Quantum Information Theory. M. R. was supported by 
DOE CSGF, Grant No. DE-FG02-97ER25308, and the 
NSF. The work at Lincoln Laboratory was sponsored by 
the AFOSR under Air Force Contract No. FA8721-05- 
C-0002. Opinions, interpretations, conclusions and rec- 
ommendations are those of the author (s) and are not 
necessarily endorsed by the U.S. Government. 



* bylander@mit.edu 

^ Present addresses: Harvard Physics Department, Cam- 
bridge, MA; 

* CIN2 ICN/CSIC, Bellaterra, Barcelona, Spain; 
§ MIT EECS Department 

[1] J. Clarke and F. K. Wilhelm, Nature 453, 1031 (2008). 
[2] R. Hanson and D. D. Awschalom, Nature (London) 453, 
1043 (2008). 

[3] Y. Nakamura, Y. A. Pashkin, and J. S. Tsai, Phys. Rev. 

Lett. 87, 246601 (2001). 
[4] D. Vion et a/., Science 296, 886 (2002). 
[5] J. M. Martinis, S. Nam, J. Aumentado, and C. Urbina, 

Phys. Rev. Lett. 89, 117901 (2002). 
[6] I. Chiorescu, Y. Nakamura, C. J. P. M. Harmans, and 

J. E. Mooij, Science 299, 1869 (2003). 
[7] M. Steffen, J. M. Martinis, and I. L. Chuang, Phys. Rev. 

B 68, 224518 (2003). 
[8] E. Lucero et a/., Phys. Rev. Lett. 100, 247001 (2008). 
[9] Y. Nakamura, Y. A. Pashkin, and J. S. Tsai, Nature 

(London) 398, 786 (1999). 
[10] T. Duty, D. Gunnarsson, K. Bladh, and P. Delsing, Phys. 

Rev. B 69, 140503(R) (2004). 



5 



[11] S. Poletto et a/., arXiv:0809.1331 (2008). 

[12] A. Izmalkov et a/., Europhys. Lett. 65, 844 (2004). 

[13] W. D. Oliver et a/., Science 310, 1653 (2005). 

[14] M. Sillanpaa, T. Lehtinen, A. Paila, Y. Makhlin, and 

P. Hakonen, Phys. Rev. Lett. 96, 187002 (2006). 
[15] C. M. Wilson et a/., Phys. Rev. Lett. 98, 257003 (2007). 
[16] D. M. Berns et a/., Nature 455, 51 (2008). 
[17] A. V. Shytov, D. A. Ivanov, and M. V. Feigel'man, Eur. 

Phys. J. B 36, 263 (2003). 
[18] S. Ashhab, J. R. Johansson, A. M. Zagoskin, and F. Nori, 



Phys. Rev. A 75, 063414 (2007). 

[19] E. M. Fortunato et al, J. Chem. Phys. 116, 7599 (2002). 

[20] S. O. Valenzuela et a/., Science 314, 1589 (2006). 

[21] J. Clarke and A. I. Braginski, eds., The SQUID Hand- 
book: Fundamentals and Technology of SQUIDs and 
SQUID Systems Vol. I (Wiley, Weinheim, 2004). 

[22] Because Ao,o, Ao,i <C Ao, 2, the population change over 
one cycle is dominated by LZ transitions at Ao,2- 



